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One of the most important predictions in magnetohydrodynamics (MHD) is that in the presence 
of a uniform magnetic field 6oe|| a transition from weak to strong wave turbulence should occur when 
going from large to small perpendicular scales. This transition is believed to be a universal property 
of several anisotropic tnrbulent systems. We present for the first time direct evidence of such a 
transition using a decaying three-dimensional direct numerical simulation of incompressible balanced 
MHD turbulence with a grid resolution of 3072^ x 256. From large to small-scales, the change of 
regime is characterized by i) a change of slope in the energy spectrum going from approximately 
—2 to —3/2; ii) an increase of the ratio between the wave and nonlinear times, with a critical ratio 
of Xc ~ 1/3; iii) a modification of the iso-contours of energy revealing a transition from a purely 
perpendicular cascade to a cascade compatible with the critical balance type phenomenology, and 
iv) an absence followed by a dramatic increase of the communication between Alfven modes. The 
changes happen at approximately the same transition scale and can be seen as manifest signatures of 
the transition from weak to strong wave turbulence. Furthermore, we observe a significant non-local 
three-wave coupling between strongly and weakly nonlinear modes resulting in an inverse transfer 
of energy from small to large-scales. 

PACS numbers: 52.30.Cv, 52.35.Bj, 47.27.Ak, 47.27.ek, 


Introduction. Waves are ubiquitous in natural sys¬ 
tems. Although waves are a basic phenomenon well un¬ 
derstood for years, the nonlinear behavior of a large en¬ 
semble of waves is still the subject of intense research 
m- When a turbulent state is developed in the pres¬ 
ence of waves one may distinguish two regimes: weak 
wave turbulence and strong wave turbulence. In the first 
case, the regime can be described analytically by a clas¬ 
sical technique based on perturbative developments [7j. 
Exact solutions of the resulting kinetic equations, cor¬ 
responding to power law spectra, can then also be de¬ 
rived [HHia- In the second case, we mainly have phe¬ 
nomenological models [14H16j : among which we find the 
critical balance (CB). Since its inception, originally in 
MHD [T71 [TB], CB has become a popular model in as¬ 
trophysics and, is believed to be also applicable to other 
systems such as electron MHD, rotating hydrodynamics 
and stratified flows USMI]. In incompressible MHD, CB 
supposes the existence of a mean magnetic guide field 
0 (which will be normalized to velocity units hereafter) 
along which propagate Alfven waves in both directions 
parallel (||) to 0. Both linear and nonlinear physics are 
affected by 0 with the development of a high degree of 
anisotropy such that energy will mainly transfer, or cas¬ 
cade, in the perpendicular (_L) direction to 0 dH 1^ . 
In such a situation, the following inequality is satished 
k± fc||. As a result of this strong anisotropy, the (lo- 
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cal) nonlinear time-scale becomes Tni ^ l/{k±b)^ whereas 
the linear Alfven wave time-scale is r^, ^ l/{k\\bQ) (for 
the derivation, see the comments on Eqs. §). The latter 
time-scale can be interpreted as the duration of a collision 
between two wave packets traveling in the opposite direc¬ 
tion at the Alfven speed 6o- The characteristic transfer 
time of energy Ttr can, as far as dimensional analysis is 
concerned, be an arbitrary function of these two times 
- an additional physical assumption is therefore neces¬ 
sary to fix the scaling. This additional assumption is 
furnished by the CB conjecture which assumes that at 
all scales in the inertial range Tni ~ Tw This physically 
means that an Alfven wave packet suffers a deformation 
of the order of the wave packet itself in one collision. Two 
properties can be derived immediately from this assump¬ 
tion: (i) the axisymmetric energy spectrum is simply of 
the Kolmogorov type, i.e. f? ~ , because the 

transfer time identifies to t„/, and (ii) a non-trivial re¬ 
lationship exists between the parallel and perpendicu¬ 
lar wave numbers in the form of /c||&o ~ k±b ~ 

[HI HD; in particular, this latter identity physically im¬ 
plies that anisotropy will increase at small-scales until 
the dissipation becomes dominant. The CB regime is 
drastically different from the weak wave turbulence one 
where, in the latter, many stochastic collisions are nec¬ 
essary to modify a wave packet significantly. In the case 
of weak turbulence, we have the inequality Tni and 

the transfer time becomes Ttr ~ /t^ [nms]. This 

transfer time can be interpreted as the time that it takes 
the cumulative perturbation (assumed to accumulate as a 
random walk) to become comparable to the amplitude of 
the wave packet itself. The resulting power law spectrum 
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in MHD corresponds to k±~'^. This result, presented 
first in a non-rigorous phenomenological way [24) . was 
then derived rigorously by a perturbative theory [5SJ US] ■ 
In particular, it is shown that no transfer (cascade) is 
expected along the parallel direction, a result stemming 
from the three-wave resonance condition |22] . The neces¬ 
sary condition for the existence of weak MHD turbulence 
is that the time ratio 
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is small (<C 1), whereas in CB it is of the order of one 
1). If we substitute the weak turbulence spectrum 
b'^Ikj_ ^ k±~'^ into Eq. Q, we see that y is an increasing 
function of k± . Therefore, there exists a critical scale be¬ 
yond which the weak turbulence cascade drives itself into 
a state which no longer satisfies the premise on which the 
theory is based. The dynamical breakdown of the weak 
turbulence description is expected to be followed by a 
saturation around one of y because of the causal impos¬ 
sibility to maintain Tw ^ ry; m- This means that for 
a sufficiently extended inertial range we should observe 
the transition from the weak turbulence regime to the CB 
one I53I13I5S]. Note, however, that CB may be refined 
by introducing the local dynamic alignment of the veloc¬ 
ity and magnetic field fluctuations which corresponds to 
a modification of the nonlinear time-scale [53]. In this 
case the power law energy spectrum is expected to be 

In this Letter, we present for the first time direct evi¬ 
dence of such a weak to strong transition, by means of a 
high resolution three-dimensional direct numerical simu¬ 
lation. 

Simulation setup. The incompressible MHD equa¬ 
tions, for our simulations, in the presence of a constant 
bo are: 


dtz^ =F bod\\z^ -h • V = —VP* -1- i^sA^z* , (2) 

where V -z^ = 0, z^ = v±b are the fluctuating Elsasser 
fields, V the plasma flow velocity, b the normalized mag¬ 
netic field (b —)■ y/fiopo b, with po a constant density 
and Po the magnetic permeability), P* the total (mag¬ 
netic plus kinetic) pressure, and a hyper-viscosity (a 
unit magnetic Prandtl number is taken). We see that the 
times Tiu and are obtained respectively from the lin¬ 
ear dispersive term and the nonlinear term on the LHS 
of Eqs. ([^, assuming a balance {z~^ ^ z~ ^ u ^ b) and 
anisotropic (kj_ ^ fey) turbulence. 

Eqs. are computed using a pseudo-spectral solver 
called TURBO [301 El] with periodic boundary condi¬ 
tions in all three directions and with 3072^ x 256 colloca¬ 
tion points (the lower resolution being in the bo direction 
where the cascade is reduced; however, the numerical box 
is not elongated and has an aspect ratio of one). The non¬ 
linear terms are partially de-aliased using a phase-shift 
method. The initial state consists of isotropic magnetic 
and velocity field fluctuations with random phases such 



FIG. 1: Top: axisymmetric spectra of the magnetic energy at 
a given = {0,1} and integrated over fey (from 2 to 128). 
Middle: time ratio y at a given fey = {2,3,4,5}. Bottom: 
integrated magnetic energy spectrum compensated by 
The dashed line corresponds to a compensated spectrum with 
k±~’^. The vertical line marks the critical scale at which the 
transition is observed. 


that the total cross-helicity, as well as the total magnetic 
and kinetic helicities, is zero (balanced and non- helical 
turbulence). The kinetic and magnetic energies are equal 
to 1/2 and localized at the largest scales of the system 
(wave numbers k G [2,4] are initially excited). We opt for 
a decaying turbulence mainly to avoid any artefact due 
to the external forcing [35]. Our analysis is systemati¬ 
cally made at a time p when the mean dissipation rate 
reaches its maximum, for which the turbulence is fully 
developed and the spectrum the most extended. Note 
that for t > t* the spectrum experiences a smooth self¬ 
similar decay which leads to a slow drift of the critical 
(transition) scale toward higher k± as the y parameter is 
proportional to the amplitude of b. We fix 1/3 = 4 x 10“^^ 
and bo = 20. Note that initially the energy of the 2D 
modes are taken to be zero in order to favor dynamics 
dominated by wave modes. With our (isotropic) initial 
conditions anisotropy will develop such that energy will 
fill the Fourier space with k± ^ fey . 

Results. We introduce the axisymmetric bi- 
dimensional magnetic energy spectrum E^{k±,ki\) 
which is linked to the magnetic energy of the system 
through the relation = ff E^{kj_, k\\)dko_dk\\. It 
is well-known that in weak MHD turbulence the 2D 
mode (/c|| = 0 ) has a singular role since it drives the 
turbulence although it is not a wave (see eg. [33]). To 
make the distinction between the contributions of the 2D 
mode and the waves we have considered the spectrum 
E^{k±, fc||) integrated from /cy = 2 to 128 (the first plane 
fcy = 1 is found to be strongly coupled with the 2D 
mode like for enslaved turbulence [33]). The result is 
shown in Fig. [^(top). At large-scale (10 < A:j_ < 100) a 
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FIG. 2: Iso-contours (in logarithmic scale) of the bi- 
dimensional magnetic energy spectrum A 

power law k^^^^ is plotted for comparison in the region cor¬ 
responding to strong wave turbnlence (see Fig. [^. 



spectrum compatible with weak turbulence is observed 
(~ k±~^). We then see that a transition seems to occur 
at a scale k± ~ 100 beyond which the spectrum becomes 
less steep. This transition happens considerably before 
the dissipative range which appears at k± ^ 600. The 
plots of the spectra for fey = 0 and 1 are also provided 
to show that there are similar, and significantly weaker 
in amplitude than the integrated spectrum. The time 
ratio Q is shown in the middle panel of Fig. for 
different (small) values of fc||. For this evaluation of x, b 

is given by 6 = ^2k\\k±E^{k±,k\\). In all cases, we see 

that x{k±,k\\) <C 1 at the largest scales as expected for 
the weak turbulence regime. The comparison with the 
spectra described above shows that a transition occurs 
when x{k±,k^\) approaches unity (> 0.1). For fey > 4 we 
find that the higher fcy , the smaller X) as expected from 
Eq. §. Note that we do not observe at small-scales 
an extended plateau where x ~ 1 which could be 
explained by the lack of resolution, the fact that we did 
not consider the local mean magnetic field and/or the 
absence of dynamic alignment in the definition of the 
nonlinear time-scale [35]. The last plot (bottom panel) 
shows the integrated spectrum compensated by k±^^^. 
The transition is visible at kj_ ~ 100 with a change in 
slope going from approximately k±~‘^ to k±~^^^; this 
happens at Xc ~ 1/3. Fig. [^displays the iso-contours 
of the bi-dimensional magnetic energy spectrum. At 
large-scale the iso-contours are strongly elongated in 
the fc_L direction meaning that the cascade is strongly 
anisotropic. At the transition scales {k± ^ 100) a drastic 
modification appears with an increase of the parallel 
transfer and therefore a stretching of the iso-contours in 
the fcy direction. Interestingly, we may find a domain 
where the edge follows approximately a power law in 
k^2/3^ This means that the energy is mainly transferred 
along an oblique direction which corresponds to CB. 
When the dissipative scales are reached the stretching 
of the iso-contours in the fcy direction increases further 


FIG. 3: Wavenumber-frequency spectrum of the magnetic 
energy E^{k±,k^\ = 5,a;). The color map is normalized to 
the maximal value of the spectrum at each fixed k±. For 
comparison we plot k±_^^^ (solid) and k± (dash-dot). The 
vertical dotted line marks the critical scale at which the tran¬ 
sition is observed and the horizontal dotted line corresponds 
to ui/{2bo) = 5 


showing a propensity toward isotropization. 

Although Fig [2 and [^ may provide a first evidence of 
a transition from weak to strong turbulence, we want to 
find other signatures. The spectrogram (wavenumber- 
frequency spectrum) of the magnetic energy provides this 
additional information. To build such a spectrogram 
one follows, in Fourier space and over a window of time 
around t*, the quantity fcy, t) = |fca;P + |fci/P + |fczP, 

at a given fcy (fcy = 5) and a given k± (from 1 to 1536). 
We then perform a time-Fourier transform of these 1536 
signals multiplied by a Hamming window and obtain 
A^(fc_L,fcy = 5,0;). The result is shown in Fig. (the 
kinetic energy is not shown but behaves similarly^. As 
we can see, at large-scales (fcj, < 60) the signal is con¬ 
centrated on a thick band localized around a;/(26o) = 5; 
thus the variable uj is closely related to fcy like the Alfven 
wave dispersion relation loa = ^||^o (the factor 1/2 in the 
frequency normalization is due to the use of the energy, 
a square of a field). That means the mode fcy = 5 com¬ 
municates only with modes directly contiguous to it, i.e. 
fcy = 4, 6 and that the system is not able to redistribute 
the energy to a wide range of frequencies, a situation ex¬ 
pected when the resonant triadic interactions dominate 
[7] . The cascade is then strongly anisotropic with a trans¬ 
fer mostly in the perpendicular direction |33j . The thick¬ 
ness of the band can be interpreted as nonlinear broad¬ 
ening due to the weak nonlinearity effects and shows that 
the resolution in the || direction is high enough to not fall 
into the discrete regime 13 [35]. From fcy ~ 60 a dras¬ 
tic change appears: suddenly the energy spreads over a 
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FIG. 4: Transfer functions T{k,p) of the total energy nor¬ 
malized to the maximum absolute value of each wavenumber 
scale. The vertical and horizontal lines mark the critical scale 
at which the transition is observed (see Fig. |^. 


wide range of frequencies. This is the most spectacular 
evidence of the emergence of a strong wave turbulence 
regime in which the || cascade is not frozen and where the 
fc|| = 5 mode becomes dynamically connected to a grow¬ 
ing number of other Alfven modes {i.e. fey ^ 5) when 
one goes to higher k±. Besides this important property, 
we note in passing that the boundary which delimits the 
region where modes are dynamically connected follows 
a power law close to which could be interpreted 

as a signature of CB since the balance condition implies 
w ^ The plot shows, however, that frequencies 

are also excited below this boundary which would cor¬ 
respond to Tw > Tni (like in the solar wind [37], but in 
apparent contradiction with a previous claim m based 
on a heuristic description and the assumption of local 
interactions). At this stage, it is important to remind 
that to define the nonlinear time-scale we have implicitly 
assumed the locality of the interactions. Then, the previ¬ 
ous observation could also be interpreted as the signature 
of non-local interactions. Note finally that in the dissi¬ 
pative range (fc_L > 600) the boundary discussed above 
seems to follow a power law close to k±^ which could be 
the signature of an isotropisation. 

The degree of locality of the perpendicular cascade 
can be investigated with the shell-to-shell energy transfer 
functions defined by: 

a*i:;“(k) = ^ [T“Jk,p) - T“(k,p)] - 2^3fc®i?“(k), (3) 

P 

dtE\]^) = ^ [Tl(k,p) - T„\(k,p)] - 2^3fc®A''(k),(4) 

P 


where [55WD) 

r#^(k,p) = ^3m{[k.Z(p)][Y(q).r(k)]}5q+p,k, (5) 

q 

is the transfer function to the mode k of field X from 
mode p of field Z, mediated by all possible triadic in¬ 
teractions with modes q of fields Y that respects the 
condition k = p + q. fJm denotes the imaginary part and 
* the complex conjugate. To study the perpendicular 
cascade, we consider concentric cylindrical shells along 
0 with constant width on a logarithmic scale which we 
define as the region < A:_l < for the 

shells numbered 4 < n < fV, where we set fco = 4 and 
TV = 31. The sum of transfer functions of the total energy 
(T = - T“ + - T\) is displayed in Fig. § While 

direct and local energy transfers dominate with transfers 
mainly concentrated around the diagonal k_\_ = p± (a re¬ 
sult found in previous decaying MHD turbulence studies 
[3311301), one observes some inverse and non-local contri¬ 
butions connecting weak and strong modes. This behav¬ 
ior is revealed by the presence of negative and positive 
energy transfer respectively in the top-left and bottom 
right part of Fig. [^ This result is new and important 
for the theory of MHD turbulence where this type of in¬ 
teraction has never been considered in the past. 

Discussion. The transition from weak to strong wave 
turbulence when passing from large to small-scales is be¬ 
lieved to be a universal property of several anisotropic 
turbulent systems with different underlying physics [211 
El- To our knowledge - and despite its importance - 
this phenomenon has so far never been observed in di¬ 
rect numerical simulations of MHD. This has left crucial 
questions unanswered. As a result of the simulation con¬ 
ducted here, we provide in this Letter a direct valida¬ 
tion of this cornerstone of anisotropic turbulence theory 
in the MHD case, and are now able to provide answers 
to some of these fundamental questions. It appears suf¬ 
ficient that the parameter xik±,k\\) crosses the critical 
value ~ 1/3 for a given fey plane to contaminate rapidly 
the others whatever their respective degree of nonlinear¬ 
ity. The spectral index and anisotropy of the total energy 
after the breakdown of the weak turbulence description 
is consistent with the establishment of CB. In addition, 
our results indicate that the transition involves blend¬ 
ing and interaction between weakly and strongly nonlin¬ 
ear modes. This unexpected behavior revealed by the 
presence of non-local and inverse energy transfers sug¬ 
gests that the transition is not simply the juxtaposition 
of weak and strong wave turbulence as was thought until 
now. This result is potentially important for other sys¬ 
tems where a transition form weak to strong wave tur¬ 
bulence is expected [7]. 
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